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ABSTRACT 



The purpose of this paper is to examine long wave propagation 
in a shallow channel of triangular cross section. Solutions for 
small scale (f = O) , large scale (-F J«? 0^ / symmetric and asymmetric 
channels are obtained. Results are shown to be consistent with 
earlier work for the fundamental mode = o}, and with edgewave 
solutions over a gently sloping bottom for the higher modes. A second 
class of waves (quasigeostrophic waves) is also obtained when the 
Coriolis effect is included. The results are compared with those for 
a rectangular channel. 
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I. INTRODUCTION 



Investigations of long wave propagation in channels with 
cross sections which don't vary along the channel began with 
Kelland [l839], who found solutions for progressive waves in a 
uniform channel of triangular cross section with sides inclined at 
45° to the horizontal. MacDonald [1894^ examined wave propaga- 
tion in a triangular channel with sides inclined at 30° to the hori- 
zontal. Most investigators have limited their work to small scale 
channels (i.e. the Coriolis effect was neglected as in the above 
cases) in which velocities normal to the axis of the channel were 
zero. An exception is Dronkers [l9643, who studied waves in a 
rotating rectangular channel, allowing cross-channel motion. 

Many rivers, estuaries, and shallow seas can be approximately 
represented by a shallow triangular channel (at least this is more 
realistic than a rectangular channel) . The shallow-water solutions 
for long waves progressing along such a channel are closely related 
to the well known solutions for edgewaves [Ursell, 1952; and Reid, 
19581 ^ Such a solution is carried out in this paper and the results 
compared to the work of Dronkers. Most of the results are obtained 
numerically. When possible, they are compared to asymptotic work. 
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II. ANALYSIS 



Long waves are studied, propagating along a channel with a 
uniform triangular cross section. The coordinate system is shown in 
Figure -I . 




Figure 1 



The maximum depth of the channel is 


ho . The slopes of the channel 




sides are s for x>o and ns 


for X < 0 (under the shallow 




water approximation s,ns<cl )• 


The depth then is 




hex') = - S X 


( x ?-o') 


(1) 


KCx') = ho + H5X 


(x < 0)^ 


(2) 



The equations to be satisfied are the shallow-water equations of motion 

-If- -f. = 0 ,3) 

= 0 , (4) 



and the continuity equation 



where u = "x - component of velocity 



V = y - component of velocity 
^ = surface height 



and y is along the channel. A wave motion in the y direction may 



be prescribed by the substitutions 

^ ^ \ <-('M g + 

^ = -7^ (x) e ' ’ ^ 



w = uCx)e 
V = V W ^ 



I 4 <r’0 



( 6 ) 



where ~tn is the wave number, and cr the frequency. Then for 



U 




(rn-fv + 



dJL.\ 

dr ) 



V = 






d-n '\ 

Jr ) 



(7) 



from equations (3) and (4). Substitution into the continuity equation 
yields 






( 8 ) 



for all r and any h(r) . 

Equation (8) will now be dealt with in the regions X > 0 and x< 0 
separately, and the results will be patched to insure wave-height and 
velocity continuity at X = 0 . 



A. EQUATION FOR X > 0 

Let I = h»/s — 'X • Then 4^X = ~ ^1, , and 

hCx) =h»-S'X = S^ . The solution region is shown in Figure 2. 
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^ =h.A 




Solution 

Region 



Figure 2 

Then (8) becomes 

CP «) 

where P = ^ _ (10) 



B. EQUATION FOR X < 0 

Let ~ y. + ho/nS . Then li6f)=-|io-+ n S;K'=n and the corres- 
ponding solution region is shown in Figure 3. 




Figure 3 



Then (8) becomes 



p' <^"n 

^ dr‘ 






where 



+ ^22if 

g 0 6' o- 



(11) 

( 12 ) 



12 



C. SOLUTION 



Using the substitutions 

■n«)= 

-n(S')^ e-^*V 



Cy > o') 
Cx > O) 



(13) 



equation (9) and (11) are put into hypergeometric form (Slater [^196(^ ); 

Jw =0 (x>0') 









(14) 



(15) 



where W~ ’X^'yn , and v/' = 3. f W • 
These equations have solutions of the form 



G — /\ Fa. C^') + (j H Cw^ C X > o') 

g' = S Fat* (w') + DH (W*) ( X < o>) 



(16) 



where 



a. 



3 



P* 

a on 



(17) 



But H, H' (u in Slater's notation) are unbounded as ^ and f 
become small and must be discarded. Then the solutions of (14) and 
(15) become 



f = /) e' ^ ^ /v (wO C X > o) 

F^'(W^) Cx< o') 
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where 



pa 



4 aW 



-^^Tf—W + C3l7~^ 



Fa'Cw') = 1 + aV + 



a'Ca' + l') ^ ,.a «'(Ta'+ OCa'-* a") . ,,J 

(i!V“‘ W W • ' ♦ ^ 



The quantities f\ and B are constants; we can set /\ = 1 as the 
results will be independent of wave amplitude in this linear theory. 

At y — 0 the patching conditions are : 

^CH)= ^CS') 0(^0 (w =W., W’= W.') (20) 

(i.e. , the wave-height is continuous), and 

U = U' At y^o (W==W., U'=Uo) (21) 



(i.e. , the cross-channel velocity component is continuous). Note 
that V/c = Tl'rnh./s and wJ = ^'»nh,/ns - W»/n are convenient 

nondimensional wave numbers. 

Using (20), (21) becomes 




Then, using (7), (20), and (23), V is also continuous at X = 0 . 
For fa, Cw’k') =■ O , (20) gives; 1) B = ^ which means (21) 
cannot be satisfied or 2) f^> (wJ) = ^ which corresponds to a 

node at X ~ O. For f^. (V/f) 0 , equations (20) and (23) give, 

eliminating B / 
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r. (w.) F.' (^) - f. (w.) ^ (-^) - F.- (^} - ^(W.)= o 




(24) 



Note that the solutions Fo^ (W») = Fa (wj) - O are contained in 
(24). Also, when the channel is symmetric (n = 1,0), (24) becomes 

Cw;')[fvrw.') - 2 4ft (w.)] = 0 ^ 

c4 W 



(25) 
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III. SMALL-SCALE CHANNELS 



For relatively short waves in narrow channels , where the 
Coriolis effect may be neglected , (10), (12), and (17) become 



P = 






= n 



and 






=r / — 3 «. 



Equation (27) can be written 

j- *'^1, 



ex. = 






= 0 - 






where << is a convenient nondimensional frequency. Similarly, 

_ X 



and 



n 



= 1 - 2 ^' 

ho _ /•. Wo 

= (i- 



Equating like terms in (28) and (30) yields the following relationship 
between cc and a' •• 

('1-2 a.') = 



Then equation (24) relates nondimensional wave numbers Wo and 
frequencies c< . 

A. SYMMETRICAL CHANNEL ( H = /. o) 

1 . Asymptotic Results 

For channels with sides of equal slope (31) gives cl-=- 
If W. < < 1 / snd a. is bounded, an asymptotic result can be 



( 26 ) 

(27) 

(28) 

(29) 

(30) 

(31) 



a. . 
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obtained. Dropping terms of OCW.^) and higher in (25) gives 



(l - - W. [ a* + a 0^ ( W. < < l') 



( 32 ) 



so that a. ~ • Thus, using (28), (32) gives 




(33) 



or 



8-< ^ Wo 



(Wo < < i) 



(34) 



Since wave celerity c = , this can be written 



= /a g 



(35) 



which is identical to Kelland's result for a symmetric triangular channel. 
Since the average depth is ho/'X / then this result also corresponds 



Other asymptotic results for V/o <<1 and a. unbounded can 
be obtained by iteration. An example is given below. 

Let Wn be the root obtained by keeping terms through 
f<K (w) ~ 0 . Then: 

- a. W, = I 



to the classical result for long waves over a flat bottom: g h . 




ck Wj — I 



a Cfk l) X I + » / - 

— rrrrr ^ 



(an'* ^ (3TF 




(xW<^ ~ etc. 
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This converges rapidly to 



- CL Wo = I. H 2 4 0 (-jc) 

and from (28) , 

~ /.*/:). ( 36 ) 

Similar techniques could give the other asymptotes suggested by the 
numerical results presented below. 

For W. >> 3 (short waves), it is expected on a physical basis 
that the waves in the two parts of the channel become independent of 
one another, and that Ursell's edgewave result 

+ O = ( 37 ) 

would hold. Although this has not been confirmed here analytically, 
it is borne out by the numerical results. 

2 . Numerical Results 

A Fortran IV computer program to solve (24) for small scale 
channels (f = o') is given in Appendix A. The results for the 
symmetric channel (y] ■=. i.o) are shown in Figures 4 and 5. In Figure 5, 
the dispersion relation for the fundamental mode = o) and the 

next five modes are plotted. 

For Wo < < 1 (Figure 4) the calculations agree well with the 
asymptotic results (34) and (36). 

For Wo >> 1 > note that 

o< ~ C + 0 (^/^- '**} 

which is the same as (37) . 
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Figure h 

Dispersion Relation: Symmetric Channel (W, << 1 ) 

(f=0,n=1 .0) 
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Figure 5 

Dispersion Relation: Symmetric Channel (f=0,n=1.0) 
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3. 



Wave Profiles 



Figure 6 is a sketch of the wave profiles for modes 
0 through 5 showing the nodal points for the symmetrical channel. 
Note that the nodes for all modes are symmetric about the center of 
the channel. The profile for the fundamental mode x o) has no 
nodes along the channel. Wave mode (rA - l) is the first asym- 
metric profile, corresponding to the solution ~ = ol- 

Wave mode has a symmetric profile with two nodal points. 

The higher modes have symmetric profiles for even mode numbers 
and asymmetric profiles for odd (fA^') , with the number of nodal points 
equal to (M). The nodal points (except 'x-o) move toward the 
channel walls with increasing Wo • This corresponds to the edge- 
wave condition for short waves (large ) . 



B. ASYMMETRICAL CHANNEL 
1 . Asymptotic Results 

For Wo < < 1 a. bounded , (24) can be written 

I - 

Since Wo<<l, then <x + a!~l. Using (31), this gives 

+ Cs-CL-t- I 

so that 

a- ^ ~ 'A • 

Let a. = ^ - € and a!~ — € * . Then (38) gives 

^ niV f ^ ) 



(38) 



(39) 
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Wave Profiles; Symmetric Channel (n=l»0) 
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which with (28) and (30) gives 

W!o — J o4 . 

This is the same as (34), and is independent of n. 

Other results for V/o << J similar to (36) are expected for 
asymmetric channels but are not examined here. 

2 . Numerical Results 

The results for n = 3.0 and n = r. o are shown in 
Figures 7 and 8. In each case, the asymptotic result (34) is confirmed. 
In fact, the entire dispersion relation for the fundamental mode (M = o') 
appears to be independent of n . The higher modes have all been 
shifted to higher frequencies in comparison with the symmetric case. 

3 . Wave Profiles 

Profiles for i.o and n =i y.O are shown in Figures 9 
and 10. Again the fundamental mode is unchanged, but the 

higher modes have been altered significantly from the symmetric case. 
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Dispersion Relation: Asymmetrical Channel (f=0,n=2.0) 
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Dxspersion Relation: Asymmetrical Channel (f=0,n=5.0) 
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Figure 10 

Wave Profiles: Asymmetric Channel (n=5*0) 
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IV. LARGE-SCALE CHANNELS 



For large channels where the Coriolis effect must be considered, 
(10), (12), and (17) give 



and 



where 




So- = 









= <r/f . 



/ 




( 40 ) 



Equating like terms in (40) , 



and 



CO- 



n ) 

r[/U'~yOf 




( 41 ) 



(42) 



where 

Now Vn 4 must satisfy both (24) and (42). 

A. SYMMETRICAL CHANNEL (n = /. o) 

1 . Asymptotic Results 

From (41), a. and a! are related by a. = a' + /or. But 
for W. << 1 and ur bounded, (42) implies that Jo.+-a'j is large, 
so that CK. -zz oJ • Then (40) gives 

-!«.= (43) 



X = 



g 



This can be written 




(44) 
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From (36) and (44) , 




ur 



(u/. < < 



Dropping the small term gives 




( 45 ) 



2 . Numerical Results 

A Fortran IV computer program to solve (24) and (42) is given 
in Appendix B. A value of )( = 0. | is chosen, which corresponds approxi- 
mately to the dimensions of the Red Sea. The results for the symmetric 
channel are shown in Figure 11. Note that for Wo < < 1 / the asymp- 
totic result (45) is confirmed. For Uf* < < 1 a new class of waves 
appears. These are analogous to the quasigeostrophic edgewaves 
obtained by Reid. Figure 12 is an expanded plot of these low frequency 
waves . 

For Wo >> 1 (short high-frequency waves), the Coriolis 

effect should be negligible. Then, fora fixed W., the frequencies <r 
for f = 0 and should be the same. Thus, 




(46) 



This relation is observed to hold for Ct/", V7o>> 1 . 



B. ASYMMETRICAL CHANNEL (r\^ l.o) 



1 . Asymptotic Results 



For Wo<<l and a. unbounded, asymptotic results similar 



to (45) are expected, but are not derived here. 
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Figure 11 

Dispersion Relatr.on: Symmetric Channel (f^0,n=1.0) 
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xon: Quasigeos trophic Waves ( f = 



2 . 



Numerical Results 



The results for n = a.o and n = 5 . o are shown in Figures 
13 and 14. Note that all modes have been shifted to higher frequencies 
in comparison with the symmetric case. Also for OJ", Wo >> 1 the 
f =• o and f SfO results give good agreement via (46) . 
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Dispersion Relation: Asymmetric Channel (fx0,n=5«0) 
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V. COMPARISON OF RESULTS FOR A RECTANGULAR CHANNEL 



AND A SYMMETRIC TRIANGULAR CHANNEL 



The dispersion relation due to Dronkers for a channel of uniform 
depth (h) and width (w) is (Appendix C) 






(n = o. i, a, . . . 



( 47 ) 



It is reasonable to compare channels of equal average depth, (Figure 15) . 




Then h = h,/a and W= , and (47) becomes 







- AlV‘ 



( 48 ) 



A. SMALL-SCALE (f = O) 

For small-scale channels (f = O) , (48) becomes 

VVo^ = ~ (tA = 0^ ^ (49) 



This agrees with (34) for the fundamental mode This dispersion 

relation (49) is compared with that for a symmetric triangular channel 
in Figure 16 . 
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Dispersion Relations: Symmetric Triangular 

vs. Rectangular (f=0) 



ro.o 
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B. LARGE-SCALE (f^O) 



For large-scale channels ii; o') , (48) becomes 

^ (a1= 0, 3, . . (50) 



This dispersion relation (50) is compared with that of a symmetric 
triangular channel in Figure 17. 
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Dispersion Relations; Symmetric Triangular 
vs. Rectangular (fxO) 
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APPENDIX A 



FORTRAN IV PROGRAM TO SOLVE PATCHING EQUATION 
FOR SMALL-SCALE CHANNELS 



IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION C(36) ,C (36) ,CP(36) , DP(36) , CCP(36) ,CDP(36),DCP(36), WORK 
*(36) ,B(36),ROOTR(36),ROOTI(36) 

1 RE AD (5, 5) A 

5 FORMAT (FI 0.5) 

IF(A.EQ. -9.50)GOTO99 
S=5.0 

AP=(S+2 . 0*A- 1 . 0)/(2 . 0*S) 

WRITE(6,10)A,AP,S 

10 FORMAT (35X, 'A= ',F10. 5, 5X,’AP= FIO . 5 , 5X, 'N = ' , F5 . 2 ,/) 

C(l) = 1.0 
C(2)=A 

DO 2 N =3,12 

C(N)=C(N-1)*(A+(N-2))/((N-1)*(N-D) 

2 CONTINUE 
D(1)=A 

D(2)=A* (a+l)/2.0 
DO 3 N=3, 12 

D(N)=D(N-1)*(A+(N-1))/((N-1)*N) 

3 CONTINUE 
CP(1) = 1.0 
CP(2)=AP/S 
DO 4 N=3,12 

CP(N)=CP(N-1)*(AP+ (N-2))/(N-l)*(N-l)*S) 

4 CONTINUE 
DP(1)=AP 

DP (2)=AP* (AP+l)/(2 . 0*S) 

DO 6 N=3, 12 

DP(N)=DP(N-1)*(AP+ (N-1))/((N-1)*N*S) 

6 CONTINUE 
DO 7 K=l, 12 
CCP(K)=0.0 
DO 7 1=1,12 
DO 7 J=l,12 
IF(I+J.NE.K)G0T07 
CCP(K-1)=CCP(K-1)+C (I)*CP(J) 

7 CONTINUE 
DO 8 K=l, 12 
CDP(K)=0.0 
DO 8 1=1,12 
DO 8 J=l,12 
IF(I+J.NE.K)G0T08 
CDP(K-1)=CDP(K-1)+C (I)*DP(J) 
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8 CONTINUE 
DO 9 K=l, 12 
DCP(K)=0.0 
DO 9 1=1, 12 
DO 9 J=l, 12 
IF(I+J.NE.K)G0T09 

DC P (K- 1) =DC P (K- 1)+D (I) *C P (J) 

9 CONTINUE 
DO 11 N=l, 11 

B (N) =C C P (N) =C D P (N) -DC P (N) 

11 CONTINUE 

CALL DPOLRT(B,WORK, 10 , ROOTR, ROOTI , lER) 
WRITE (6,25) (ROOTR (I) , ROOTI (I) , 1= 1 , 1 0) 

25 FORMAT(//,4X,2D25 . 16,/) 

WRITE(6,30)IER 
30 FORMAT (5 OX, 12) 

GOTOl 
99 STOP 
END 
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APPENDIX B 



FORTRAN IV PROGRAM TO SOLVE PATCHING EQUATION 
FOR LARGE-SCALE CHANNELS 



IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION C(36) , D (36) , CP(36) , DP(36), CCP(36), CDP(36), DCP(36), WORK 
*(36) ,B(36) ,ROOTR(36) ,ROOTI(36) 

CF=1.0D-04 
H=1000.0 
G=9.8 
SLP=0. 1 
S=5.0 
CAPA=. 10 

1 READ (5, 5) W 

5 FORMAT (FI 0.5) 

IF(W.GT. 14.5)G0T099 
WRITE (6, 10) W 

10 FORMAT(50X,'W=',F10.5,/) 

OMEGA=C.50 
12 OMEGA=OMEGA-.001 

IF(OMEGA.DE, 0 . 0)GOTOl 

A= . 5 * ( 1 . 0/OMEGA+ 1.0- (CAPA* (OMEGA* OMEGA- 1.0)) /W) 

AP= . 5* (1 . 0-1 . 0/OMEGA-(CAPA* (OMEGA* OMEGA- 1 . 0)/(W*S))) 

C(l) = 1.0 
C(2)=A 
DO 2 N=3,12 

C(N)=C (N-l)*(A+(N-2))/((N-l)*(N-D) 

2 CONTINUE 
D(1)=A 

D(2)=A*(A+l)/2.0 
DO 3 N=3,12 

D(N)=D(N-1)*(A+(N-1))/((N-1)*N) 

3 CONTINUE 
CP(1) = 1.0 
CP(2)=AP/S 
DO 4 N=3,12 

CP (N) =CP (N - 1) * (AP+ (N -2) )/( (N- 1) * (N- 1) *S) 

4 CONTINUE 
DP(1)=AP 

DP(2)=AP* (AP+l)/(2 . 0*S) 

DO 6 N=3 , 12 

DP(N)=DP(N-1)*(AP+ (N-1))/((N-1)*N*S) 

6 CONTINUE 
DO 7 K=l, 12 
CCP(K)=0.0 
DO 7 1=1,12 
DO 7 J=l, 12 
IF(I+J.NE.K)G0T07 
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CCP(K-1)=CCP(K-1)+C (I)*CP(J) 

7 CONTINUE 
DO 8 K=l, 12 
CDP(K)=0.0 
DO 8 1=1, 12 
DO 8 J=l, 12 
IF(I+J.NE.K)G0T08 
CDP(K-1)=CDP(K-1)+C (I)*DP(J) 

8 CONTINUE 
DO 9 K=1 , 12 
DCP(K)=0.0 
DO 9 1=1, 12 
DO 9 J=l, 12 
IF(I+J).NE.K)G0T09 
DCP(K-1)=DCP(K-1)+D(I)*CP(J) 

9 CONTINUE 
DO 11 N=l,ll 

B (N) =CC P (N) -C DP (N) -DC P (N) 

11 CONTINUE 

FCN=B(1)+B(2)*W+B(3)*W*W+B(4I*W**3+B(5)*W**44-B(6)*W**5+B(7)*W* 
*+B(8)*W**7l-B(9)*W**8+B(10)*W**9+B(ll)*W**10 
WRITE (6 , 1 5) A , AP , OMEGA , W , FC N 
15 FORMAT(2X,4F15.7,D25.16) 

G0T012 
99 STOP 
END 
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APPENDIX C 



DISPERSION RELATION FOR DRONKERS' RECTANGULAR CHANNEL 



A. FORMULATION 

The equations to be satisfied are: 






~ "F V + q -C ^ j. a ^ £- _ 

K ( AL ^ 0 



- i y + '’■t) 



ol'r< 



Substitution of 

( u( , V, j=) = ( V. r?) e’ 

into (1) and (2) gives 

I <r u - fv + g-^ = 0 

i.o-v+fu+ 

L<rn.O. 

Then , for cr is^ -f *■ , 

LI 

V= 

and the ordinary differential equation for is 



A 



j X 



t + ^ = 



where 






with solution of the form 



= a. AAyr\ fix b cxt<u "X , 



( 1 ) 

( 2 ) 



(3) 



(4) 

(5) 

(6) 

(7) 
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B. BOUNDARY CONDITIONS 

The coordinate system is shown in Figure Cl. The flux normal 
to the channel walls must vanish at the walls. 

////// ////////////// y y / y 






x= o 






"v \ \ \ V \ \ 



Figure Cl 



Then (4) gives 



-2nf2L + ^ 



ctI-X 



Substitution of (7) into (8) yields 



(y^ oy) . 



f 



^ flW- = 0 



SO that. 



or 



^ IV = 0 

= 0,1 ,a , • .•) 



f\^ = V' 



Then (6) gives the dispersion relation 



'yri^ = 









/Vjrr-^ 

w"" • 



( 8 ) 



(9) 



(10) 



( 11 ) 
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